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I. INTRODUCTION 



The calculation of the recoil corrections of order (Za) 6 (m/M)m to the hydrogen energy 
levels has a long history Refs. Jl|-§[. After initial disagreements consensus was achieved in 
Ref. 0, where one and the same result was obtained in two apparently different frameworks. 
The first, more traditional approach, used earlier in Refs. [@-^J], starts with an effective Dirac 
equation in the external field. Corrections to the Dirac energy levels are calculated with the 
help of a systematic diagrammatic procedure. The other logically independent calculational 
framework, also used in Ref. Jj]], starts with an exact expression for all recoil corrections 
of the first order in the mass ratio of the light and heavy particles m/M. This remarkable 
expression, which is exact in Za, was first discovered by M. A. Braun ||, and rederived 
later in different ways in a number of papers Refs. PJlO|7|. 



The agreement on the (Za) 6 (m/M)m contribution achieved in JF| seemed to put an 
end to all problems connected with this correction. However, it was claimed in a recent 
work |TTJ, that the result of [[!]] is in error. The discrepancy between the results of Refs. 
|7|JTT|| is confusing since the calculation in |ll|] is performed in the same framework as the 
one employed in |7|, namely it is based on a particularly nice form of the Braun formula 



obtained by the author earlier |10 



AE rec = --^Re J ^<n\(p-t(u))G(E + u)(p-f>(u))\n>, (1) 

where summation over all intermediate states is understood, G(E + uS) is the Coulomb- 
Green function in the Coulomb gauge, which in the momentum space has the form 

D( W , k) = -4wZa(* - ^) 2 } = -AnZa 2 "* (2) 

and 

at = TV • (3) 

Note that T){cu,k) is nothing more than the transverse photon propagator with the source 
at the proton position, and integration over the exchanged photon momentum k is implicit in 
the expression above. Below we will explicitly perform multiplication in the matrix element 
in Eq. ([!]). Respective contributions to the energy levels will be called Coulomb (corresponds 
to pp), magnetic (corresponds to pD and Dp), and seagull (corresponds to DD). 

It is the aim of this paper to resolve the above noted discrepancy on the recoil correction 
of order (Za) 6 (m/M)m to the IS" energy level, and also to obtain this correction for the 
S- levels with arbitrary principal quantum number (it was earlier calculated only for n = 1, 2 



II. TWO APPROACHES TO THE BRAUN FORMULA 

Calculation of the recoil contribution of order {Za) 6 generated by the Braun formula 
was performed in [|7| in a most straightforward way since separation of the high- and low- 
frequency contributions was made in the framework of the e-method developed by one of the 
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authors earlier [0. Hence, not only were contributions of order (Za) 6 (m/M)m obtained in 
Ref. 0, but also linear in m/M parts of recoil corrections of orders (Za) A and (Za) 5 (ref. 

were reproduced for the 15-state. Note that the Braun formula, despite its obvious 
advantages, in its present form sums only contributions linear in the mass ratio. Hence, old 
methods are more adequate for obtaining the proper mass dependence of the contributions 
of orders (Za) 4 and (Za) 5 , which were worked out in Ref. |]. Calculations in Ref. @ turned 
out to be rather lengthy and tedious just because all corrections of previous orders in Za 
were reproduced. 

The most significant feature of the recoil corrections of order (Za) 6 , which made the 
whole approach of Ref. |jTT| possible, is connected with the absence of logarithmic recoil 
corrections of this order, as was proved in ||. Unlike [0], the calculations in Jll| are orga- 



nized in such a way that one explicitly makes approximations inadequate for calculation of 
the contributions of the previous orders in Za, significantly simplifying calculation of the 
correction of order (Za) 6 . Due to absence of the logarithmic contributions of order (Za) 6 , 
infrared divergences connected with the crude approximations unadequate for calculation of 
the contributions of the previous orders would be powerlike and can be safely thrown away. 
Next, absence of logarithmic corrections of order (Za) 6 means that it is not necessary to 
worry too much about matching the low- and high-frequency (long- and short-distance in 
terms of Ref. [llj) contributions, since each region will produce only nonlogarithmic contri- 
butions and correction terms would be suppressed as powers of the separation parameter. 
We would like to emphasize once more that this approach would be doomed if the loga- 
rithmic divergences were present, since in such a case one could not hope to calculate an 
additive constant to the log, since the exact value of the integration cutoff would not be 
known. 

We are going to perform below calculation of the recoil contribution of order (Za) 6 in 
the framework of Ref. |[L1|| , and to discover the source of discrepancy between the results of 
Ref. [0] and Ref. In order to really implement such program we need to have a regular 
method to qualify all terms which will be thrown away. To this end we will use a slight 
generalization of the ordinary approach to calculation of the leading order contribution to 
the Lamb shift. 

It may be proved that all corrections of order (Za) 6 (m/M)m are generated by the 
exchange of photons with momenta larger than m(Za) 2 , so we will consider below only this 
integration region. In the spirit of the common approach to the Lamb shift calculations we 
will split the integration region over the exchanged photon momenta (and when necessary 
over frequencies) with the help of an auxiliary parameter a which satisfies the conditions 

mZa <C a Cm, (4) 

and we will call the photons with momenta smaller than a low-frequency (or long- 
distance) photons, and the photons with momenta larger than a will be called high-frequency 
(or short-distance) photons. Considering low-frequency photons we may expand over the 
ratio k/m since for such photons k/m < a/m <C lQ. On the other hand, for the high- 



1 Note that the apparent linear divergences in this region of the form a/m are really parametrically 
small. 
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frequency photons mZa/k < mZa/cr <C 1, and we may expand over this parameter. Note 
that for momenta of order a both expansions are valid simultaneously, and, hence, we may 
match the expansions and get rid of the auxiliary parameter a. However, the problem under 
consideration is in a sense even simpler than calculation of the leading order contribution to 
the Lamb shift, and due to absence of the logarithmic contributions of order (Za) 6 (m/M)m, 
precise matching of the high- and low-frequency contributions is unnecessary. Below we will 
consider calculation only of the low-frequency (mZa < k < a) contribution to the energy 
shift, since for the high frequency contribution the results of Ref. ]7) and Ref. |~T| nicely 
coincide. 



III. MAIN RECOIL CONTRIBUTION 

With the help of the Braun formula one may easily obtain an expression for the leading 
recoil correction which is linear in the mass ratio and which includes all terms of order (Za) 4 
and lower (see Ref. To this end we rewrite the Coulomb contribution in Eq.([l|) in the 
form 

AEcoui = ^ < n \p 2 \ n > -jj < n|pA~p]|n > (5) 
= AE cl + AE c2 . 

We also extract the nonretarded Breit part from the magnetic contribution in Eq.(|l] ) 

^^'inagn AEg r -\- AE ma g n r , (6) 



where 



and 



AEBr = < n|pf)(0 ' k) + f)(0 ' k ^ n >' (7) 



A£ m(W = ~J ^ < n\[V,p]G(E + cu)-b(cu,k) (8) 



-f)(u,k)G(E + u)[V,p]\n > ' 



u + iO' 



where V is the Coulomb potential (V = —Za/r). 

Now it is not difficult to check with the help of the virial relations (see, e.g., Ref. ||14]| ), 
that the sum of the main part of the Coulomb term and of the Breit contribution acquires 
a very nice form 

AE cl +AE Br = m2 ~^\ (9) 

where E is the value of the energy given by the Dirac equation. As we will see below, 
all other recoil contributions to the energy level start at least with the term of order (Za) 5 , 
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and, hence, the formula above correctly describes all contributions of order (Za) 4 and lower. 
However, this formula describes only contributions linear in the mass ratio. A more precise 
expression which takes into account corrections of higher order in m/M, was obtained in 
Ref. 0. 

It is easy to see that the expression in Eq.(|9[) also contains the correction of order (Za) 6 , 
which for the nS-states has the form 

13 1 1 AZafm . , 

AE GY = (- + - j + n^Ti" 1 ' 10 

This contribution was originally obtained in Ref. [|TJ. 

The remaining part of the Coulomb contribution has the form 

AE c2 = --^ < n|pA_p|n > . (11) 

Let us check that this term leads to corrections of higher order than (Za) 6 when the 
intermediate momenta are of the atomic scale. We want to exploit the large (of order 2m) 
value of the energy gap between positive and negative states in comparison with the typical 
energy splittings (of order m(Za) 2 ) in the positive energy spectrum. First, let us note that 

< n|[p,V]A_[p,1/]|n >=< n\[p,H - E]A_[p,H -E]\n > (12) 
= — < n|p \m >< m\(E n — E m ) 2 ~p\n > . 

However, (E n — E m ) 2 > 4m 2 (1 — ca 2 ), and, hence, 

| < n|[p, V r ]A_[p, V]\n> \ = \ < n|p^ \m >< m\(E n - E m ) 2 p\n > \ (13) 



> | < n|pA_p|n > |4m 2 (l — ca 2 ). 

Then 

| < n| P A_p|n > | < t _ | <7i|[p,V]A-[p,V]|n> 1- (14) 

We know that at the atomic scale the Coulomb potential is of order (Za) 2 , the momentum 
operators are of order Za, and, hence, we explicitly have the factor (Za) 6 . Note that this 
approach would not work if we had a projector on the positive energy states. In such a 
case the energy differences would be of order (Za) 2 themselves and we would not get any 
suppression, since the factors (Za) 2 would cancel in the numerator and denominator. 

Returning to our case, it is easy to realize that the projector on the negative energy 
states leads to additional suppression in the nonrelativistic limit, and, hence, the term 
under consideration does not produce any contribution of order (Za) 6 at the atomic scale. 

There is complete agreement between the results of Ref. [|7| and Ref. [11| for the correc- 
tions discussed in this section. 
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IV. SEAGULL CONTRIBUTION 



Following Ref. |TTJ let us again start with the Braun expression Eq.(|T]) for the seagull 
contribution and perform the integration by closing the contour each time around one of the 
transverse photon poles 

AE S = ~ f^-.< n|D(w, k)G(E + cj)D(cu, k)\n> . (15) 
M J 2m 

Substituting the pole representation for the Coulomb Green function we obtain in accor- 
dance with Ref. [pj 



a p ( Za ) 2 ^ | 47TQfe> rv ^ \m >< m\ . E m -E ^ 

s ~ 2M 1 k> L V (E-k'-E m )(E-k-E m ) [ + k> + k> 1 ' 



^ (E + k> - E m )(E + k - E m ) [ k' + k )l k |n> ' 

Let us consider positive- and negative-energy parts of this expression separately. 

We may expand the positive energy part in (E — E m )/k and (E — E m )/k, taking into 
account that in the low-frequency integration region mZa < k < a. In the first order of 
this expansion we get 

a 7-i+ ( Za ) 2 „ ,47rafc/ 4vrOfc 

= W <n| ^ A+ ^ |n> ' 

Calculation of this contribution will be considered below. Let us turn to the negative- 
energy contribution. Energy differences are large for the negative energy contribution (\E — 
E m \ « 2m(l — ca 2 )), so we expand the negative energy term in kj (E — E m ) 

V \m><m\ E m — E 

_ (E + k' — E m )(E + k — E m ) k' + k' K ' 

Elm >< ml r 1 (k + A;') 2 n 

J L 1_ 



^ £-£ m 'fc + A;' 2(E-E. 



3 J 



In accordance with Ref. |TTJ the terms linear in k/2m cancel, and the negative energy 



contribution acquires the form 

(Za) 2 Airct k > , 1 (k + k 1 ) 2 n 47re* fcl , , 

AE~ = - - < n\ —A I - ■ - - n> (19) 

4mM{l + ca 2 ) 1 k' k + k' 2[2m(l + ca 2 )} 3 k 1 ' 1 ' 

It may be shown (compare below consideration of the negative energy contribution in 
the case of the one transverse exchange) that the first term produces the contribution of 
order {Za) 5 while the second is of order {Za) 1 . Only terms linear in k, k' are capable of 
producing contributions of order {Za) 6 , but these terms cancel each other, as we have just 
seen. 
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Let us now return to the positive energy contribution. The idea of Ref. [[LI]] is to consider 
matrix elements and to calculate them in the nonrelativistic approximation, which produces 
the leading low-frequency contribution. All matrix elements under consideration have com- 
mon structure. In general they are the products of matrix elements of 7-matrices in the 
momentum space. Each such matrix element in the nonrelativistic limit may easily be re- 
duced to an explicit function of momenta and cr-matrices, then transformed into coordinate 
space and calculated between Coulomb- Schro dinger wave functions. 

We have performed an explicit calculation along these lines and obtained in complete 
accord with Ref. [pH 

. 1 (Za) 2 ,1 1 31 2 + 2<tZ, 

AEt = ) n { < n 2p— p + — — A n>. 20 

4m 2 M 1 r 2 ^ r 4 2r 4 1 v ; 

This expression is singular at the origin. This singularity produces linear and logarithmic 
ultraviolet divergences in momentum space as well as a constant contribution, and hence the 
contribution under consideration cannot be calculated unambiguously in the general case. It 
is necessary to realize at this stage that the initial expression for the seagull contribution in 
Eq. (|TSp was defined unambiguously. Even separation of the integration region with the help 
of the auxiliary parameter a could not lead to an ultraviolet divergence in the low-frequency 
region since all momentum integrations are cut off from above by o and should generate 
not power divergent but power suppressed terms. It is clear that the apparent divergence 
is connected with our inaccurate calculation of the singularity at large momenta or small 
distances. Hence, we have to return to the initial momentum space expression for the positive 
energy seagull contribution and perform all calculations directly in the momentum space. 
The result of such a calculation may be later interpreted as an unambiguous prescription 
for the proper regularization of the coordinate space operators for the .S-states. 

Note, that for the non-S* states, wave functions vanish at the origin, the operators above 
are well defined on such wave functions, and lead to unambiguous results. Of course, any 
regularization at small distances will not influence the value of the non-S* matrix elements 
of the operator in Eq. (|20|) , and will not influence the agreement between the P- level energy 
shift calculated in Ref. ||11|| , and the same shift obtained earlier in another framework in Ref. 



A. Accurate Calculation with Momentum Space Cutoff 

Direct calculation of the positive energy seagull contribution Eq.(|16D in momentum space 
leads to the following expression for the S'-state contribution 

2 „ ^3^1 A3 m a3ij j3u 0^2 



(Za) 2 r d 3 p' d 3 p d 3 k' d 3 k 



71 



AE - == 5f J ^Wf^fWf^ p ' ~ p " k " k 'fe < 21 > 



Ar , (k'k)(p'k')(pk) k'k 1/N . 

V(P')[-P P + % k 2 ~ —MP) = AS -1 + AE s2 + AV" 



The first two terms in the integrand do not rise too rapidly with k and k', and we may 
unambiguously calculate them using the Fourier transforms discussed above. For the first 
term we have 
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A F - ( Za ^ f^ r d3 P' d "P d3k ' ^(-p'+p+k+k'H 87T " (02) 

sl ~~ m 2 M J (2tt) 3 (2tt) 3 (2tt) 3 M 3 ' k' 2 k 2 1 ' 



m?M J (2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 

^(p')pW(p) = / ^^^ iP( " p,+p) )>^)pWb). 

The remaining integration over p' and p simply returns us to the coordinate space wave 
functions, and we may rewrite the expression above in the operator notation^ 

= bib < n|p ^ > • (23) 

This contribution exactly reproduces the nonsingular operator obtained in the previous 
chapter. 

Next we calculate the second contribution in the same manner as above 

32 ~ m 2 M J J (2vr) 3 (2tt) 3 (2tt) 3 (2tt) 3 W ' W W 

(24) 

(Za) 2 f , 3 f d 3 p' d 3 p ir{ , } >\P'jP m (Z ^iur \ 



Now we use the formula 

(rp')(rp) = -[r x p'][r x p] + r 2 (p'p), (25) 

and omit the terms with the vector product since we are considering only ^-states now. 
Then we obtain 

AE s2 = 0. (26) 

Next we have to calculate the third contribution, which corresponds to the 1/r 4 term 
in the naive result above in Eq. (PU|) . This time we cannot use Fourier transformations over 
exchanged momenta for calculation of this integral, since this leads to a singular expression 
in coordinate space. So we first perform the safe Fourier transformations over the wave 
function momenta, and then directly evaluate the exchanged momenta integrals, taking into 
account that they are cut from above by cr <C m, 



2 One has to take into account that the apparent sign of the expression below changes, since the 
momenta in the exponent have opposite signs. 
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(27) 



In order to preserve the transparency of the presentation we will perform the calculation 
only for n — 1 here. The general case of arbitrary principal quantum number will be 
considered at the end of the paper. We substitute explicit expressions for the lS-wave 
functions in the formula above, and do the coordinate-space integral 

^ - -^p/SS^/^^ (28) 

QAnHZaf . .,.„/■ d 3 k' d 3 k k'k 
-7l^(0)| 2 1 



m 2 M " rv 71 J (27r) 3 (27r) 3 fc /2 P[(k + k') 2 + (27) 



212 ' 



where 7 = mZa. 

Symmetrical integrals over the exchanged momenta are cut from above by the parameter 
a. However, first integration, say over k', is convergent at high momenta and the cutoff may 
be safely ignored 

. _,,_4 IQn(Za) 2 . ,, f d?k f°° ,,, f 1 , k'kx , , 

^ -^t^' /n^I (29) 



tW0)| 2 t^^[ — r^ x -H = -^rf7^(0)| 2 dfc[ 



m 2 M " rv " 7 (2tt) 3 P l fc 2 7 J m 2 M " rv " Jo L fc 2 7 

-Tl^(0)| 2 [27T In ■ 



m 2 M 27 7 J 

Nonlogarithmic term of order (Za) 5 in this expression is additionally suppressed by the 
small ratio a/m, and may be safely ignored. Thus, we see that the properly regularized 
operator 1/r 4 in the seagull diagram does not generate a constant contribution. The loga- 
rithmic divergence above should cancel with the respective contribution of the one-transverse 
(magnetic) diagram. 

V. MAGNETIC CONTRIBUTION 

This time we start with the Braun expression for the one transverse photon in Eq.(|l]) 

AE magn = ±-Re I < n\pG(E + u)f>(u, k) + k)G(E + u)p\n > (30) 
and first calculate the contour integral! 3 ] 



3 Note that the overall minus sign is connected with the respective sign in the definition of the 
transverse propagator. 
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a n -Za i \ m >< m \ \ - I rri >< ml n 47ra: fcl . . 

AE maan = < n p > > — : — \n>+h.c. 31 

As we are again calculating the low-frequency corrections to the Breit potential let us 
expand the positive energy term in (E m — E)/k 

&E+ agn (32) 

Za „ 1 E m -E (E m — E) 2 Ana k 

= — < ra p > \m >< m h • • • — ; — pi > +n.c. 

2M IH Y fc fc fc 

The first term in this expansion may be written in the form 

A^+ a5nl = < n| P A + -^|ra > +/I.C. (33) 

Za 47ro:fc Za 47ro:fc 

< n, p — —5 — |?t. > H — — < ra|pA_ — — — n > +n.c. 



2M ^ k 2 1 2M P 



A£ Br + AE + 



magnl— ) 

and it is now evident that the first (Breit) term here coincides with that part of transverse 
exchange which cancels with the respective term in the Coulomb contribution. 
Remaining positive-energy contributions are given by the expression 

1 Za , m ,r E m — E (E m — E) 2 Attolu. 

&E+ agnr = -— < n|p£ \m >< HI" V + 1 fc3 + • • -Inr I" > + h - c - ( 34 ) 



A -I- A 4- 



A. Positive Energy Contribution 



In accordance with Ref . Jl l] one may check that the term Ai?^ a „ n2 does not lead to the 



contributions of order (Za) . We have 

Za * -, 

2M <n|P V'"" " lv "' m " y fc 3 



A£+ affn2 = ^ <n|p^|m><m|(£ m -£)-^|n>+/i.c. (35) 



(Za) 2 ,47rfc' A Ana kl 

< n\——A + —-^\n > +h.c. 



2M 1 k' 2 ^ k 3 

The simplest way to estimate this matrix element is to make a Fourier transformation. 
Then we need an infrared divergent Fourier transform of 1/k 3 . All momentum integrals in 
the low-frequency region are cut off from below by m(Za) 2 , and it is easy to check that the 
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leading term in the infrared divergent Fourier transform generates a logarithmic divergent 
contribution of order (Za) 5 in accordance with Ref. The next terms vanish with the 
infrared cutoff and cannot produce contributions of order (Za) 6 . 

Let us turn now to the term AE^ agn3 . Naive calculation in the coordinate space in 



accordance with the result in Ref. [Ill leads to the result 



&E+ agn3 = < n\ P J2(E m - Ef\m >< m\^£±\n > +h.c. (36) 



+ 



Za) 2 . 1 71 2 aL 

< n\2p— p - — - —\n > 



Am 2 M r 2 2r 4 r 

This expression contains only operators which are nonsingular at the origin for ^-states. 
Hence, they are well defined, and there is no need for a careful momentum space considera- 
tion in this case. 



B. Negative Energy Contribution 

There are two negative-energy contributions connected with the magnetic term, one in 
Eq.@Tp, and the other in Eq.(||). 
Let us consider first 

. Za x , \m >< ml Anct k 

^E~ agn = — < n\ P Y, E _ Em + k k I- > +h.c (37) 



We have checked, in accordance with Ref. Ill, that this term leads at most to contribu- 



tions of order (Za) 6 , and, hence, is of no interest. 

We still have to calculate one more negative energy contribution, contained in Eq. (|33|) 

&E+ agnl _ = — < n| P A_-^|n > +h.c. (38) 

= (Zal Mc*4*Hak) 
8m 2 M 1 k' 2 k 2 1 

Naive calculation with the help of the Fourier transformation leads, in accordance with 
Ref. to the expression 



Ai Wa- " Am 2 M < Hl r > ' [6J) 

However, this expression, as in the case of the seagull contribution, contains singular 
operators at the origin, and does not have unambiguous meaning for the S-states. A more 
careful calculation, which explicitly takes into account a momentum space cutoff a, is needed. 

First we transform the negative energy contribution in Eq.(^8|) to the form 

AE+ ma9nl - = < n\[p, V}A.^\n > +h.c. (40) 
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Next we substitute the negative energy projection operator in the nonrelativistic approx- 
imation A_(p) m 1/2 — {cap + (3m) /2m and use the trivial identity 



[p, V}A_ = A_[p, V] - [A_, [p, V]] = A_[p, V] + [p, V]]. (41) 

Note that the first term on the right hand side vanishes applied to the ket-vector, and 
the negative energy contribution reduces in the nonrelativistic approximation to 

AEL ^- = -2^u\ (2^ < ^\P k \p,V]^\n > . (42) 



Then we use 



< rc(r)|p k = -27 < n(r)\y, (43) 
[p,V] = -i(Za)^, 

and obtain 

» (^") 2 /" ^ /" rt N2 (r fc r)47re ikr 

A£ -»- = bar / csf / d r « r » V— • < 44) 

As in the case of the singular seagull contribution we will perform the calculation for 
n — 1 first, postponing consideration of the general case to the next chapter. We substitute 
explicit expressions for the wave functions and obtain 



ab -»>- - w*<°>' 2 iw&L dx{i - x2) r <45) 



MZaf a 2 r d 3 fc 47r r (47 2 + fc 2 )arctan 



-7IVW 



2 J 



2m 2 M ' Irv " 7 (2vr) 3 A; 2L 7 fc 3 fc 



47r(Za) 2 (4tt) 2 _ 2l/ _^ |2 [° JLr ( 4 7 2 + fc 2 ) arctanj^ _ 9 

2 J 



2m 2 M (2VT) 37 ^ (0)l J ^ lk s k 

Again, as in the case of the seagull contribution, this term may be understood as a proper 
regularization of the naive singular in the coordinate space operator from Eq.([3~9l). 
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VI. CALCULATIONS FOR ARBITRARY PRINCIPAL QUANTUM NUMBER 



The total low-frequency contribution for the lS'-state is given by the sum of the results 
in Eq.(0), Eq.fl23|), Eq.(||), Eq.(||) and Eq.fllF 



AE, 



low—freq 



'is) 



(46) 



and coincides with the result obtained earlier for the low-frequency contribution in Ref. 
|7|]. We see that the seagull and magnetic contributions partially cancel each other. This 
reflects cancellation of the 1/r 4 terms in the language of Ref. |TT| . However, the contribution 
(— 1) survives. This contribution is connected with the 5-function term in Ref. [O], and the 



error in Ref. [I I ] is due to an improper regularization of this contribution. Note that from 
the point of view of the coordinate representation after the Fourier transformation is done 
the proper regularization is highly nontrivial. One could never obtain this contribution with 
a naive ad hoc regularization in coordinate space. 

The result in Eq.(|46|) is valid only for the lS'-state. We are going to generalize it to an 
arbitrary principal quantum number. 



A. Seagull Contribution for Arbitrary nS'-Level 

The general expression for the wave function of an nS-level has the form 

T 3 i Tl — 1 

Mr) = (^) 5 e"-[l ir + ...]. (47) 

Let us introduce (3 = j/n. Almost all calculations above for n = 1 immediately turn 



into calculations for arbitrary n after substitution j — > (3 [ 16[ . The wave function has the 
form 

Mr) = (— )5e-^[l - (n - l)(3r + ...} = ^ n (0)e^ r [l - (n - l)(3r + ...]. (48) 

7T 

Quadratic and higher order terms in r in the postexponential factor in the wave function 
do not produce any contribution to the energy level connected with the singular operator 
in the naive expression in Eq.(|20|), and we will ignore them below. The only difference 
between the general case and the case of n = 1 is connected with the linear term in the 
postexponential factor. Let us find out how it changes the result for the seagull contribution. 
First, let us write down the singular seagull contribution induced by the purely exponential 
part of the wave function for arbitrary n in the form 

AEl »' = -^\ m p fpLjL 4 J^0l f^We-m (49) 

m 2 M J (27r) 3 (27r) 3 k' 2 k 2 J 



m 2 M 
where 



{ -^\m\ 2 ei /r ' 
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l /r 4 _ „, O 



< /r = 2 ^ ln ^ ( 5 °) 
The linear terms in the wave functions lead to an additional contribution 

A ^ = -^i^(°)i 2 / 1^ (C^ 1 ™^ / ) re -^[_ 2(ri _ 1)/3r] (51) 



m 2 M in " 1 " dp s m 2 M' rv 71 v n ' 2(3 

and the total seagull contribution to the energy shift is equal to 

AEl& = + A^ rr = -^l\ m \^r + (n - 1)/?^^] (52) 

W,. ( .MM 2 



|^(0)| 2 [27T7ln— -27r(n- 



m 2 M' rv ' 2/5 



B. Magnetic Contribution for Arbitrary nS-Level 

As in the case of the seagull contribution the only difference of the general case from the 
case of n — 1 is connected with the linear term in the postexponential factor in the wave 
function. The purely exponential part of the wave function leads to the following singular 
magnetic contribution for arbitrary n 

The new contribution induced by the linear term in the wave function has the form 

Za r d 3 k 47re ikr 
&E+ agnl _ >cor = -^2^ J (27F [ " (n " 1)P]{< n|rPfc[Pj V] H^ ]n > (54) 

47re ikr 

+ < n\p k [p, V] k2 r\n >}. 

Next we write 

rp k = p k r - [p k ,r], (55) 
and using the commutation relation 



[Pk,r] = -iy, (56) 



obtain 
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. j_ Za r d 3 k r . , ^. r , r k r . 47re* kr , , , 

A^,.^ = -—/^-jjK,,- l)fl { < [P, l'] — In > (57) 

47re* kr 

+ < w|Pfc[p,V] fc2 2r\n >}. 
Za f d 3 k r k T/1 4vre ikr 

= "wm! (2^ (n " 1){< -K-^7 )[P ' > 

47re ikr 

-(n - < n\p h \p,V]— p— 2r|n >} 

= ^ l^(0)| 2 (- " l)[4 a9 nx- +/3 2 ^(^)] 

Then the total singular magnetic contribution is equal to 

^Kna g m~,tot = &E+ agnl _ + AE+ agnl _, cor = i^|^( )| 2 [2vr/31n ^ - tt/3 (58) 

+27r/3(n - l)ln^ - 3(n - 1)tt/3] = ^^(0)| 2 [27r 7 In ^ - tt/? - 3(n - 1)tt/3]. 



VII. TOTAL RECOIL CORRECTION 

The total low-frequency contribution of order (Za) 6 (m/M)m for arbitrary nS-state is 
given by the sum of the terms in Eq.fllOD, Eq.(^3|), Eq.(|52|), Eq.(|36|) and Eq. (|58|) 

,13 1 1 AZaf m (Za) 2 , 1 . 

AEiow-freq = (- + - - - + —) — -m + — < „|p_p|„ > (59) 

[ Z(t) #(0)| 2 [27r 7 ln - 2n(n - 1)(3 - 2a] - < n|2p^p|n > 



m 2 M irv 2/3 v ' J 4m 2 M 1 r 



+ ^ IV;(0)|2[27r7ln 2^ _ ^ _ 3(n _ 1)7r/?] 
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,13 1 1 AZaf m (Zaf m 

— ( I I p L 777, _ J: L 777 

v 8 8n n 2 2n 3> n 3 M n 3 M 

Note that the last term connected with the naive singular operators in the coordinate 
space turned out to be state-independent. 

To obtain the total recoil correction of order (Za) 6 (m/M)m it is also necessary to cal- 
culate the high-frequency (or short-distance) contribution to the energy shift. The simplest 
way is to use again the Braun formula Eq.(^), but this time in the Feynman gauge. This 
calculation is quite straightforward if one again uses the auxiliary parameter a introduced 
above in order to qualify would be infrared divergences. Such a calculation was performed 
explicitly in ref. ]TT]] and led to the result 



AE high . freg = (4 In 2 - ^)^L^m, (60) 
in complete agreement with Ref. |7]]. 

Then total correction of order (Za) 6 (m/M)m to the energy levels is given by the sum 
of the results in Eq.floTID and Eq.flBOD 

13 1 1 AZaf m , , lAZafm , . 

AE tot = (- + - j + ^P^-Trm + 41n2 - -)^-J-— m . 61 

8 8n n z In 6 n 6 M 2 n d M 

For n = 1, 2 this result nicely coincides with the one obtained in Ref. 0. 

In conclusion let us emphasize that discrepancies between the different results for the 
correction of order (Za) e (m/M) to the energy levels of the hydrogenlike ions are resolved 
and the correction of this order is now firmly established. 
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